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Differential coronal rotation using radio images at 17 GHz 
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ABSTRACT 

In the present work, we perform time-series analysis on the latitude bins of the solar 
full disk (SFD) images of Nobeyama Radioheliograph (NoRH) at 17 GHz. The flux 
modulation method traces the passage of radio features over the solar disc and the 
autocorrelation analysis of the time-series data of SFD images (one per day) for the 
period 1999-2001 gives the rotation period as a function of latitude extending from 
60°S to 60°N. The results show that the solar corona rotates less differentially than 
the photosphere and chromosphere, i.e., it has smaller gradient in the rotation rate. 
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1 INTRODUCTION 

In the recent studies of the solar rotation, the tracer methods 
were extensively used for the determination. The features 
like sunspots, plages, filaments, faculae, bright points, super 
granules, coro nal holes, giant cells, etc. were found to rotat e 
with the Sun (|Howardill99ll . Il996l : IShivraman et al.lll993l ). 
By tracing the passage of any of these features over the solar 
disc, one can estimate the rotation as well as its differential 
nature with respect to height and latitude. In general, the 
equatorial region rotates faster than the polar-region. 

The solar coronal rotation is a comparatively less un- 
derstood phenomenon because the coronal features are less 
prominent in both the time duration and spatial extent. The 
observations are affected by the fact that the solar corona is 
optically thin across a wide range of observed radio frequen- 
cies. Recent studies of the solar radio emission provided a 
lot o f inform ation concerning the solar corona (|Vats et al.l 
ll998l : lMehtal [2005V The radio observations have the advan- 
tage of permitting one to attain resolution for the heights in 
the corona as the propagation of radio waves distinctly de- 
pends on the plasma density in th e corona. By usin g radio 
flux data at different frequencies (IVats et aLlbooJ l. it has 
been shown earlier that the coronal rotation depends on the 
heights in the corona. 

Several groups have reported from time to time that the 
solar corona rotates differentially, with respect to latitude, 
in the same way as the photosphere and the chromosphere. 
The ana lysis of X-ray bright points ( XBPs ) in SOHO/EIT 
images (Braisa et al.l l2001al . l2002al . 12004 iKarachik et al.l 
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|2006| ). and recent a nalysis of XBPs in Hinode/XRT and 
Yb/ifco/i/SXT images l|Karivappall2008h also asserts the coro- 
nal diffe rential rotation. Ho weve r, in contrast to these ref - 
erences, IWeber et all l| 19991 ) and lWeber fc Sturrock I (|2002D 
studied the rotation rate of corona by using Yohkoh/SXT 
data and showed that the soft X-ray corona does not exhibit 
any significant differential rotation, i.e. the corona rotates 
more rigidly than the chromosphere and the photosphere. 

Recently, by making use of XBPs observed with the 
Yohkoh/SXT, iHaral (|2009l ) found that the rate of differen- 
tial rotation changes with a parameter At which is linked 
with the lifetime of XBPs. The long-lived XBPs follows the 
rotation rate of the photospheric magnetic field, whereas, 
the short-lived XBPs approaches to that evaluated from the 
photospheric Doppler measurements. 

IZaatri et all (|2009l ) compared differential rotation of 
subphotospheric layers derived from GONG-I— I- doppler- 
grams with the small bright coronal structure (SBCS) ob- 
served through SOHO/FIT for the period correspond to the 
decling phase of solar cycle 23. They reported that the SBCS 
rotate faster than the upper subphotospheric layer (3Mm) 
by about 0.5 deg/day at the equator. The latitude gradient 
of the rotation rate of the SBCS and the subphotospheric 
layers were found to be very similar. 

The comparison of the data of the solar magnetic el- 
ements reveal that the differential rotation of the com- 
pact magnetic elements with negative and positive polar- 
ities have similar beh avior for the solar cycles 20 and 21 
(|japaridze et al.ll2009l ). It is also established that some vari- 
ations in the rotation rate of compact magnetic elements 
were present at the time of polarity reversal of the Sun. 

The differential rotation of the solar corona is, there- 
fore, still an open subject and further studies are required 
to establish the coronal differential rotation. Here we report 
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Figure 1. Typical autocorrelogram of the year 1999, used in the present work for the estimation of synodic rotation period. 



the study of the coronal rotation using Nobeyama Radiohe- 
liograph (NoRH). 



2 OBSERVATIONS AND METHODOLOGY 

The Nobeyama Radioheliographs at 17 GHz consist of an 
array of eighty four antennas aligned in the T-shaped con- 
figuration. This is dedicated only for the solar observations. 
High speed processing of signals from all antennas generates 
full disc radio images of the Sun, at the maximum rate of 20 
images per second. The Sun is observed for 8 hours per day. 
The radioheliograph provides unprecedented uniform sam- 
pling of SFD images with high spatial and te mporal resolu - 
tions of 10 arcsec and 50 msec, respectively l|Shinzdll995h . 
The observations are possible throughout the year without 
any interruption by weather conditions. The radio images 
seem to discriminate low-contrast or faint structures on the 
solar disk e.g. dark filaments, bright points, coronal holes, 
and other large-scale structures. 

Here, we analyze the solar full disk images obtained 
from Nobeyama Radioheliograph (NoRH) at 17 GHz. Each 
image is of 512 x 512 pixel in size. The principal technique 
of our investigation is to calculate the autocorrelation co- 
efficient of the intensity variation for a set of latitude bins 
selected on the NoRH images. For this we utilize one NoRH 
full disk image per day on which latitude bins, or strips, are 
chosen. The dimensions of these rectangular bins are just 
two pixels in width and the length is such as to incorporate 
the entire pixels on the solar disc image at each latitude. 
The strips are selected at every 10° latitude 

For each latitude, a time series of the radio intensity is 
generated from the average of the selected pixels values in 
the strip. To determine rotation period, we obtained auto- 
correlation coefficient using the standard subroutines of IDL 



Table 1. The coefficient A & B obtained from Nobeyama 
RadioheUograph (NoRH) at 17 GHz for the years 1999- 
2000. The last column gives annual sunspot numbers. 



Years 


Coefficient 


Annual 




A±Ea 


B±Eb 


Sunspot 








Numbers 


1999 


14.826 ± 0.200 


-1.528 ± 0.474 


093.3 


2000 


15.083 ± 0.239 


-2.892 ± 0.567 


119.6 


2001 


14.543 ± 0.097 


-1.968 ± 0.229 


111.0 


Mean 


14.817 ± 0.060 


-2.129 ± 0.143 


107.9 



software. The period estimation is a straightforward identifi- 
cation of any dominating periodically varying signal present 
in the coronal radio emission. 



3 DATA ANALYSIS 

For the time-series of each latitude, autocorrelation upto a 
lag of 145 days are calculated and plotted in the form of au- 
tocorrelogram. The typical autocorrelogram for a data set 
for the year 1999 at 30° north latitude is shown in Figure 1. A 
systematic periodic modulation can be seen from the auto- 
correlogram as a modulation due to the solar rotation. The 
autocorrelation coefficient at some of the latitudes are quite 
high and the curve is quite smooth showing many damped 
cyclic oscillations (Figure 1). The amplitude of oscillatory 
features reduce perhaps due to the temporal variability in 
the solar rotation. This is because of the long term persis- 
tence of the solar features on the solar disc at that latitude. 
The short term features at higher latitudes affect adversely 
the smoothness, numbers of cyclic variation as well as the 
magnitude of autocorrelation coefficient. Nevertheless, the 
first peak of each autocorrelogram is quite clear and could 
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Figure 2. The coronal sidereal rotation period as a function of 
the latitude for the year 1999 
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Figure 3. As in Figure 2 for the year 2000 



be used for the estimation of the synodic rotation period 
upto ±60° latitude. The average radio flux at higher lati- 
tudes (> 60°), in both the hemispheres, does not show pe- 
riodic nature and hence is not used for obtaining rotational 
periods. Thus from NoRH images we are able to get rota- 
tional periods from 60° S to 60°N only. The radio emission 
at higher latitudes is either uniform or totally random as 
a function of time. Employing this method, we determined 
first the synodic rotation period and then from the synodic 
the sidereal rotation period, at every 10° latitude in the range 
of ±60° each year during 1999-2001. 

To compare the sidereal rotation rate Q{tp) of this work 
with other observations and their models of the solar rota- 
tion, we flt the data using the standard polynomial expan- 
sion in its most commonly used form, as given below 



n{i>) =A + Bsm'ip 



(1) 



where f2(V') is the sidereal rotation rate at heliographic 
latitude ip. The parameter A represents the equatorial ro- 
tation rate and B represents differential rate largely at low 
latitudes. Since we have not used the data of the higher 
latitudes > 60° on both the sides of the equator, hence, it 
is sufficient to use onl y the first t wo terms of the standard 
polynomial expansion (|Beckll 19991 ). 

The rotation rate is related to the rotation period as 
given below. 



360° 



(2) 



where T{ip) is in day and in degree/day. 



4 COMPARISON 

The estimated sidereal rotation period has been plotted as a 
function of the latitude for the years 1999-2001, separately. 
Figures 2-4 represent the differential rotation profiles of 
each year. In spite of scatter in the data points, the solar 
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Figure 4. As in Figure 2 for the year 2001 



differential rotation is clearly evident. We obtained the co- 
efficient A and B of the solar coronal rotation for each year 
separately. These are tabulated in Table 1. The mean of 
these has been compared with the results obtained through 
the analysis of Hmode/XRT XBPs, Yohkoh /SXT XBPs 
(lKarivappall2008l ) and SOHO/EIT CBPs data (|Braisa et all 
|2004 ). 

The sidereal equatorial rate (coefficient A, in Table 2) 
in the case of NoRH is found to be 4.4 % higher than the 
Hinode/XRT XBPs (14.2 deg/day) and 15.8 % lower than 
the Yohkoh/SXT XBPs data (17.6 deg/day), on average. 
Coefficient A of the present work is found to be compara- 
ble with the sunspot data (14.6 deg/day) and SOHO/EIT 
CBPs data (14.7 deg/day). The differential rate represented 
by coefficient B is f ound to be s ignific antly lower then the 
values obtained by iKarivappal (120081 ) using Hinode/XRT 
XBPs data (49.4 %) & Yohkoh/SXT XBPs data (53.1 
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Table 2. The coefficients A &i B obtained from the present work and other recent pubhshed results. 



Data Source 


Tracers 


Coefficients 


References 


Period 






A±Ea 


B±Eb 






NoRH 


flux modulation 


14.817±0.060 


-2.129±0.143 


present work 


1999-2001 


Hinodc/XRT 


X-ray bright points 


14.192±0.170 


-4.211±0.775 


Karivappa (2008) 


Jan, Mar & Apr 2007 


Yohkoh/SXT 


Coronal bright points 


17.597±0.398 


-4.542±1.136 


Karivappa f2008) 


1992-2001 


SOHO/EIT 


Coronal bright points 


14.677±0.033 


-3.100±0.140 


Braisa et al. ('2004) 


Jun 1998-May 1999 


Greenwich 


Sunspots 


14.551±0.006 


-2.870±0.060 


Balthasar et al. { 1986) 


1874-1976 
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Figure 5. The coefficients A, B and annualy averaged sunspot 
numbers for 1999-2000. 



%). lBraisa ||2004) usinR SOHO/'EIT CBPs data (31.3 

%) and lBalthasar et al.l l| 19861 ) using photospheric sunspots 
data (25.8 %). It means the latitude gradient of the coro- 
nal rotation rate is much lower at 17 GHz radio frequency in 
comparison to that obtained from the observations of X-Ray 

and EUV frequencies. 

From the analysis of radio flux at 2.8 GHz. lMehtal (|2005l ) 
could not find any systematic relationship between coronal 
rotation period and phases of solar cycle. Our results show 
comparatively faster equatorial rotation in the year 2000, 
than those of 1999 and 2001. Thus, it may have some rela- 
tion with the maxima of the solar activity in the 23d cycle. 
The sunspot numbers in 1999 and 2001 are relatively lower 
than 2000. The estimated coefficients A and B from the 
Nobeyama radioheliogram show a temporal variability. The 
temporal variations of A and B and sunspot numbers are 
shown in Figure 5. 



5 CONCLUSION 

Our results for coefficient A are more in agreement 
with the resul ts derived from the sunspot region at the 
photospheric l| Howard et al. Il984l : iBalthasar et al.l 1 19861 : 
IPulkkinen fc TuominenI Il998l) and at the chrom ospheric 



level ijBraisa et al.l l2004l : iKarachik e'tal\ l2006h . Since 



Yo/ifco/i/SXT observed the soft X-ray emission of the corona 



and was not a coronograph, the portion of corona that was 
imaged through this SXT range from the photosphere to a 
height of approximately 0.5 Rq (David McKanzie, Private 
communication 2009). Whereas, the radio emissions at 17 
GHz originate approximately at the altitude of ^ 1.2 x 10* 
km above the p hotosphere. The height is adopted from 
(|Vats et al.ll200l[). wherein the enha nced electron density 
model by Aschwanden fc Bend l| 19951 ) was used to estimate 
the average altitude. Thus 17 GHz seems to originate almost 
at the interface of the chromosphere and the corona, hence, 
the equatorial rotation rate matches more with the photo- 
spheric or the chromospheric data rather than the coronal 
data. The analysis of daily solar full disc radio images at 17 
GHz for the years 1999 to 2001 shows that the coronal equa- 
torial rotation rate is higher than those obtaine d by XBPs 
and quite lower than those obtained by CBPs (|Karivappal 

Hoos). 

The coefficient B shows that the corona does rotate dif- 
ferentially as in the photosphere and chromosphere; but the 
gradient is lower (< 25.8% and <31.3%, respe ctively). This 
is in agree ment with the results obt ained by IWeber et al.l 
(|l999h and lWeber fc Sturrockl (|2002h . The difference in the 
rotation profiles of soft X-ray emission and 17 GHz could 
possibly be due to their different regions of origin in the so- 
lar atmosphere. The coefficients A and B are also compared 
with the annual sunspot numbers (see Table 1). This gives 
evidence of the dependence of coronal rotation on the phases 
of solar cycle. 

Figure 5 shows that the coronal equatorial rotation (^4) 
is in phase with the annual sunspot number, whereas its lat- 
itude dependent differential gradient is anti-phase to the an- 
nual sunspot number. A long term study is indeed required 
to ascertain this correlation more clearly. 
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